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Abstract 

We use a microscopic Slater-Koster tight-binding model with short-range exchange and atomic 
spin-orbit interactions that reaUstically captures generic features of ferromagnetic metal nanoparti- 
cles to address the mesoscopic physics of magnctocrystallinc anisotropy and hysteresis in nanopar- 
ticle quasiparticle excitation spectra. Our analysis is based on qualitative arguments supported 
by self-consistent Hartree-Fock calculations for nanoparticles containing up to 260 atoms. Calcu- 
lations of the total energy as a function of magnetization direction demonstrate that the magnetic 
anisotropy per atom fluctuates by several percents when the number of electrons in the particle 
changes by one, even for the largest particles wc consider. Contributions of individual orbitals to 
the magnetic anisotropy are characterized by a broad distribution with a mean more than two or- 
ders of magnitude smaller than its variance and with no detectable correlations between anisotropy 
contribution and quasiparticle energy. We find that the discrete quasiparticle excitation spectrum 
of a nanoparticle displays a complex non-monotonic dependence on an external magnetic field, with 
abrupt jumps when the magnetization direction is reversed by the field, explaining recent spectro- 
scopic studies of magnetic nanoparticles. Our results suggests the existence of a broad cross-over 
from a weak spin-orbit coupling to a strong spin-orbit coupling regime, occurring over the range 
from approximately 200- to 1000-atom nanoparticles. 
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I. INTRODUCTION 



Interest in the properties of magnetic nanoparticles has grown recently, partly because of 
advances in synthesis and measurement techniques and partly because of potential applica- 
tions for high storage-density magnetic media and spin electronics. Ferromagnetic nanopar- 
ticles with diameters of a few nanometers containing of order of 1000 or fewer atoms can 
now be reliably fabricated and studied with a variety of different methods|]l|, 0, Small 
monodomain magnetic particles have traditionally^, |^ been described using classical mi- 
cromagnetic theory, in which the total energy is expressed as a function of the magnetization 
orientation. Shape and magnetocrystalline magnetic anisotropy leads to a dependence of 
energy on orientation, to barriers that separate minima that occur at easy magnetization 
orientations!^, and to hysteretic discontinuous changes in orientation as a function of the 
strength of an external magnetic field. When the size of a magnetic particle is only a few 
nanometers, the discrete nature of its quantum energy spectrum can be directly observable 
at low temperatures and starts to affect the magnetic properties of the particle. 

A milestone in the experimental study of individual ferromagnetic metal nanoparticles 
was achieved recently by Deshmukh, Gueron, Ralph et.a/.[|T], |7|, using single-electron tun- 
neling spectroscopy. By exploiting the Coulomb blockade effect, these experimentalists were 
able to resolve individual quantum states in the discrete many-body excitation spectrum 
of single ferromagnetic metal nanograins with sizes in the range from 1 to 4 nm. As in 
bulk ferromagnetic metals, the low-energy excitations of the nanoparticles that were probed 
in these experiments involve both particle-hole excitations of the electronic quasiparticles 
and quantized collective excitations of the magnetization-orientation field that appears in 
the classical micromagnetic theory. In an initial attempt to achieve an understanding of 
the novel physics evident in the external field dependence of the experimental excitation 
spectra, two of us recently § analyzed a simple quantum model with long range exchange 
interactions. We were able to demonstrate explicitly that the low-energy excitations of a fer- 
romagnetic metal nanoparticle are specified by the occupation numbers of its quasiparticle 
orbitals, as in a Fermi liquid, and by the global orientation of the aligned spins of all single- 
occupied (and therefore spin-polarized) orbitals. This model is not, however, able to account 
realistically for the influence of spin-orbit interactions, which play the essential role in con- 
trolling the complex hysteretic behavior seen in these experiments |^. Kleff, von Delft, 
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Deshmukh and Ralph p, ^ have proposed that the single-electron tunneling spectroscopy 
experiments can be explained by accounting for non-equilibrium spin accumulation, and by 
assuming that the magnetic anisotropy energy of a small magnetic particle has surprisingly 
large fluctuations as a function of the number of electrons on the particle. This assumption 
leads to a non-trivial magnetic field dependence of tunneling resonances that resembles ex- 
perimental behavior. More critically, these authors point out that if non-equilibrium spin 
and quasi-particle excitations both occur |]lOl, the low energy spectra are characterized by 
many closely spaced resonances, consistent with experiment 0] . In Ref. we presented a 
possible approach toward achieving a unified and consistent quantum description of both 
collective and quasiparticle physics in magnetic nanoparticles. 

The attempts put forward in Refs. H, ||, |10| to develop a quantum description of ferromag- 
netic metal nanoparticles do not address the microscopic origin of the magnetic anisotropy 
and the fluctuations of this quantity that are a necessary consequence of spin-orbit inter- 
actions in mesoscopic systems. This article addresses these issues and, more in general, 
investigates the changes in magnetic properties of small metallic nanoparticles that occur 
because of the finite spacing between quasiparticle levels near the Fermi energy. Our con- 
clusions follow from qualitative arguments based on perturbation theory expressions for the 
influence of spin-orbit interactions on quasiparticle energy levels, and on numerical studies 
of a simplified model that we believe is sufficiently realistic to describe generic aspects of 
the interrelated mesoscopic physics of quasiparticle energy-levels and magnetic anisotropy 
energy in ferromagnetic metal nanoparticles. We are particularly interested in the variation 
of quasiparticle energies with fields on the scale of the coercive field, which is trivial in the 
absence of spin-orbit interactions, but entirely non-trivial in their presence. The model we 
study is based on a Slater-Koster tight binding Hamiltonian and a mean-field treatment of 
exchange interactions. Our aim is to understand the dependence of magnetic anisotropy 
and quasiparticle energy-level-spacing statistics on particle size and shape, external mag- 
netic field, and, with single-electron-transistor systems in mind, also on electron number. 
We find that, because of the absence of strong correlations between angular momentum op- 
erator matrix elements and orbital energy differences, the size of spin-orbit induced energy 
shifts in nanoparticles and bulk perfect crystals are similar. On the basis of estimates for 
energy shifts and for avoided crossing gaps, we predict that a crossover from weak to strong 
quasiparticle spin-orbit scattering will occur over the range from approximately 200-atom 



4 



to approximately 1000-atom nanoparticles. We find that for small particles the contribution 
of individual quasiparticle orbitals to the magnetocrystalline anisotropy energy has a wide 
distribution, characterized by a variance comparable to the spin-orbit-scattering lifetime 
broadening energy, Tit^q, and a mean that is smaller by more than two orders of magni- 
tude. Surprisingly, we find no measurable correlation between the contributions to magnetic 
anisotropy from quasiparticles that are close in energy in this limit. As a result of the statis- 
tical properties of the quasiparticle magnetic anisotropies, the total magnetic anisotropy per 
atom fluctuates by several percents when the number of electrons in the nanoparticle changes 
by one, even for particles containing 260 atoms. Finally, in agreement with experiment 0, 
we find that the quasiparticle excitation spectrum exhibits a complex non-monotonic behav- 
ior as a function of an external magnetic field, with abrupt jumps when the magnetization 
orientation of a nanoparticle changes discontinuously in response to the field. Our analysis 
provides insight into mesoscopic fluctuations of the micromagnetic energy functional which 
appears in the classical theory of a magnetic nanoparticle. 

Our paper is organized in the following way. In Section H we introduce the model and 
describe the formalism. In Section |TT1| we analyze the qualitative change in quasiparticle 
energy-level statistics induced by spin-orbit interactions, and discuss the connection between 
the quasiparticle properties and the magnetic anisotropy of a ferromagnetic nanoparticle. 
Fluctuations of the magnetic anisotropy as a function of spin-splitting field, atom number, 
and electron number are investigated in Section |IV|. Magnetic hysteresis and the exter- 
nal field dependence of quasi-particle energy levels are discussed in Section 0. Finally, in 
Section we summarize our findings and present our conclusions. 



II. THE MODEL 

We model the nano-particle as a cluster of Ma atoms located on the sites of a truncated 
crystal. The numerical results we present here are for a cobalt cluster whose truncated f.c.c. 
crystal is circumscribed by a hemisphere whose equator lies in the x|/-plane of the f.c.c. 
crystal fS^ . 
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A. Tight-binding Hamiltonian and Slater-Koster parameters 



The model we use is intended to qualitatively capture the physics of a transition-metal, 
itinerant-electron ferromagnet. We use a s-p-d tight-binding model for the quasiparticle 
orbitals, with 18 orbitals per atom including the spin-degree of freedom. Nine orbitals per 
Co atom are occupied for in neutral nanoparticles. The full Hamiltonian is, 



'H — 7-^band + 'Hcxch + '^SO + '^Zee • (l) 

Here Hband is a one-body term describing the orbital motion of the electrons. In second 
quantization, it has the form 

"^band ~ 5Z S ^ lli,fi2,s^l,tJ.i,s^j,tJ.2,s (2) 

i,j Mi,M2 

where and c are Fermion creation and annihilation operators for single-particle states 
labeled by i, fi, s. The indices i,j are atomic site labels, and f'^ couples up to second nearest- 
neighbors. The indices /ii,/i2 label the nine distinct atomic orbitals (one 4s, three 4p and 
five 3d). The spin degrees of freedom, labeled by the index s, double the number of orbitals 
at each site. It is useful for us to vary the spin-quantization axis, which is specified by a unit 
vector f2(G, $) where B and $ are the usual angular coordinates defined with respect to the 



f.c.c. crystal axes. The parameters t^-J s Slater-Koster parameters ||12|| obtained after 



performing a Lowdin symmetric orthogonalization procedure ||13|] on the set of Slater-Koster 
parameters for non-orthogonal atomic orbitals of bulk spin-unpolarized Co||T4|. 

The electron-electron interaction term in Eq. (|l|) is simplified by introducing explicitly 
||38| only the ferromagnetic exchange interaction, Hexch, between the electrons spins of d- 
orbitals on the same atomic site. These interactions are largely responsible for magnetic 
order in transition metal ferromagnets: 



'^exch — —'^Udd X! ■ Sd^i , (3) 

i 

where 

fi(^d fi^d s,s' 

The parameter Udd in Eq- (0) determines the strength of the exchange interaction and is set 
equal to 1 eV. 
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This value of Udd gives rise in our finite clusters to an average magnetic moment per 



atom of the order of 2/iB, which is larger than the bulk value in Co|0, in agreement with 
other calculations [|I5| and experimental result s|l6|, p!7| for Co clusters. This value of Udd is 
also approximately consistent with a mean-field relationship, derived below, between band 
spin-splitting and magnetization, using bulk values]!^ for these two quantities. In Eq. (p, 
r is a vector whose components r"^, a = x,y,z are the three Pauli matrices. 

The third term in Eq. (||), Tiso, is a one-body operator, essentially atomic in character. 



representing the spin-orbit interaction. It can be written as [18 



n 



so 



^dJ2 J2 • 5|/i',s') 



l,lM,S"t,ll' ,s' 



(5) 



The atomic matrix elements {fi, s\L- S\fi' , s') = {i, fi, s\L- S\i, ^' , s') can be been calculated 
explicitly as a function of the direction of the magnetization ri[^. The energy scale ^d, which 
characterizes the coupling between spin and orbital degrees of freedom, varies in the range 
from 50 to 100 meV in bulk 3d transition metal ferromagnets|pO|| . For our calculations we 
have used ^d = 82 meV, taken from Ref. pT| . 

Spin-orbit coupling gives rise to a dependence of the total energy of a ferromag- 
net on the direction of its spontaneous magnetization, an effect known as magneto- 



crystalline anisotropy||39||; with spin-orbit coupling the magnetization is partially orbital 
in character [^, ^ and is sensitive to orbital anisotropy due to crystal- field interactions on 
an atomic site, due to the spatial arrangement of the atoms neighboring that site to which 
electrons can hop, and due to the overall shape of the full nanoparticle that becomes avail- 
able to an electron after many hops We will see below that the hemispherical shape 
of the nanoparticles we have studied plays the dominant role in determining their magnetic 
anisotropy energy. 

The last term in Eq. (|l]) is a local one-body operator, representing the Zeeman coupling 
of the orbital and spin degrees of freedom to an external magnetic field -ffext^ 



'Hzee= -^bXI Y1 {f^^s\{L + gsS\fl',s')) ■ He^tclf,^sCi^^',s' 
i fj,,fj,' ,s,s' 

= ~f^B ^ext • { I] (/^, s\L\fl', S)) cl^^sCi,f,',s + ^ 4,f,,s^s,s'Ci,^,s'} ■ (6) 

The extreme sensitivity of magnetization orientation to external magnetic field is a com- 
bined effect of the collective behavior of many electrons in a ferromagnetic nanoparticle and 
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the smallness of the magnetic anisotropy energy relative to the overall energy gain associated 
with ferromagnetic order. The most delicate physics of a ferromagnetic nanoparticle, and in 
our view the most interesting, is that associated with magnetization direction reorientation 
by weak external magnetic fields. Thus the interplay between the Zeeman term and the 
magnetic anisotropy produced by spin-orbit interactions is at the heart of the physics we 
intend to address. 



B. Mean-field approximation 

We seek a ferromagnetic solution to the mean-field equations for this model, decoupling 
the quartic term in the exchange interaction using the ansatz 

Sd,i = {Sd,i) + SSd,i , (7) 

ignoring terms that are second order in SS^^i , and determining ground state expectation 
values, (...), self-consistently. This standard procedure leads to a Hamiltonian that can be 
diagonalized numerically and to a self-consistency condition that can be solved iteratively to 
determine the mean-field order parameters. For the present calculation we have simplified 
this procedure further, by averaging the spin-splitting exchange mean field, 

h, = hn = -^2{Sd,) (8) 

over all sites. Our motivation for doing so is to simplify the magnetic anisotropy energy 
landscape discussed below, forcing all spins to change their orientations coherently. We rec- 



ognize that complicated non-collinear spin configurations p^, |25| commonly occur in mag- 
netic nanoparticles, and that under the action of an external field small groups of atoms can 
change their orientation relative to other parts of the nanoparticle. Complex magnetiza- 
tion reorientation processes are obvious in addition spectroscopy experiments!]^. By forcing 
atoms to change their magnetic orientations coherently, we are restricting our attention to 
relatively large nanoparticles (A/^ > 50) in which most atoms have parallel spins and to 
nanoparticles with simple bistable hysteretic behavior in which the physics we address will 
be easier to study. The beautiful series of detailed SQUID magnetometry experiments by 
Wernsdorfer and colleagues 0, ^ demonstrate that nanoparticles can be prepared that do 



have simple coherent magnetization reversal properties. 
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By using a simplified model Hamiltonian we are able to deal with larger nanoparticle 
systems than would be possible with a first principles calculation [jl5|, ^ since we 
are interested only in generic aspects of the ferromagnetic nanoparticle physics there is little 
to gain from the additional realism that could be achieved by performing self-consistent 
spin-density functional calculations. 

Given the averaged spin-splitting field h, the mean-field Hamiltonian is now a single-body 
operator 

h-h ^ ^ 

T-CuFih) = TYband + ^SO + ^Zee + 7^7— (gsUB)^ Na - Qsl^Bh ■ 2Sd,i ■ (9) 

The self-consistent spin-spitting field is also the field, denoted by h* below, at which the total 
ground-state energy function E{h) = (Hmf(^)) is minimized. Notice that gsfJ^Bh* has the 
dimension of an energy. In fact, in absence of spin-orbit interaction, 2gsfiBh* can be identified 
with the band spin-splitting A = ep^ — e^^, where ep^ and ep^ are the majority and minority 
spin quasi-particle Fermi energies. |^ Diagonalization of the quadratic Hamiltonian yields 
a set of quasi-particle energies {en,s}, n = 1,2, . . ., which in the ground state are filled up to 
a Fermi energy determined by the number of valence electrons in the nanoparticle. Strictly 
speaking, because of the spin-orbit interaction, the spin character of the corresponding 
eigenstates \'ipn,s) is not well defined. 

For small nanoparticles, most eigenstates have predominantly spin-up or spin-down char- 
acter and we sometimes use this property to assign spin labels s = t for spins along the 
order direction and s = | for reversed spins. 

III. SPIN-ORBIT INTERACTIONS AND THE MAGNETIC ANISOTROPY EN- 
ERGY 

The magnetic anisotropy energy of small ferromagnetic particles § has two fundamentally 
distinct origins, long-range magnetic-dipole interactions which cause a dependence on overall 
sample shape and short-range exchange interactions that, because of atomic-like spin-orbit 
interactions, are sensitive to all aspects of the electron hopping network including bulk 
crystal symmetry, facet orientations, and also overall sample shape. We concentrate here 
on spin-orbit-induced magnetocrystalline anisotropy which gives rise to the most interesting 
physics in ferromagnetic nanoparticles. When magnetostatic shape anisotropy is important, 

9 



it can be added as a separate contribution. We begin our discussion with some qualitative 
estimates of the effect of spin-orbit interactions that are based on perturbation theory|^, |29| . 



A. Perturbation theory considerations 

In bulk 3d transition metal ferromagnets, spin-orbit interactions are relatively weak. 



Their coupling strength is less than 10% of the d-band width Wd in bulk materials pO 
allowing the energy shifts they produce to be estimated perturbatively. Because of angular 
momentum quenching in the absence of external fields, the expectation value of H^q is zero, 
even in case of ferromagnets [pO|] . 

The quasi-particle energy shift due to spin-orbit interactions is given by second-order 
perturbation theory as 



|2 



_ _ VSd; I \rr?i,s' l-^l rn,s/ 's',s\ , s 

eso = (^n,s - e„,s - 2^ _ .0 ' \^^) 

s' n,s m,s' 

where \4'n,s) e° are respectively the single-particle eigenstates and energies in the 
absence of spin-orbit interaction. In small particles, the importance of the spin-orbit inter- 
actions can be assessed by comparing the spin-orbit energy shift eso with the single-particle 
mean-level spacing 6. 

In an infinite periodic solid only states at the same k are coupled and these are separated 
energetically by an energy comparable to the bandwidth Wd- In a nanoparticle a given state 
will be coupled to many other orbitals, but the coupling matrix elements will be reduced in 
accord with the following sum rule: 

E • v,.r = ■ ^lo ~ 4. (11) 

s' 



The estimate for the right hand side of Eq. ^ is based on the atomic character of the 
angular momentum in our model, and uses that L ■ L ~ [5*6 + 3*2 + l*0]/9, with 
the estimate of the typical L ■ L expectation value representing an average over d, s and p 
orbitals. It follows that, unless there are important correlations between angular momentum 
matrix elements and quasiparticle energy differences, the typical shift in energy caused by 
the spin-orbit interaction is eso ~ both bulk crystals and in nanoparticles. 
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which is in the range between 1 to 10 meV. For example in Co, using = 82meV and 
Wd ~ 5eV, this rough estimate gives for the magnitude of the spin-orbit energy shift eso 
~ 1.3meV. The sign of the shift might be expected to be sensitive to the spacing of nearby 
quasiparticle orbitals. The anisotropy energy, that is, the dependence of the total band 
energy on the magnetization orientation, is given to a good approximation by a partial 
canceling sum of spin-orbit induced energy shift dependences on magnetic orientation. 

In the approximation that the exchange field is orbital independent, majority spin and 
minority spin orbitals are identical and differ only in their occupation numbers. In this 
approximation there is no contribution to the anisotropy energy from doubly occupied or- 
bitals. Because of the cancellations, the anisotropy energy per atom is much smaller than 
eso- For example, the zero-temperature anisotropy energy per atom in bulk is 60/ieV for 



h.c.p. Co, and ~ IfieV for b.c.c. Fe and f.c.c. Ni|^0|. In a finite or disordered system, there 
will always be perturbative coupling to quasiparticle states close in energy in Eq. |10|, but 



the matrix elements, which satisfy the sum rule of Eq. |Tl|, will be distributed among many 
states and typical energy shifts in nanoparticles should generally be comparable to those in 
bulk perfect crystals. Typically net anisotropy energies per atom in small magnetic particles 
are larger than the bulk because of the loss of symmetry at the surface. 

An important quantity used to characterize the strength of spin-orbit interactions in bulk 
systems and large nanoparticles is the the spin-orbit scattering time, rsoPT], 

In the weak coupling regime, it is given by the Fermi's golden rule 

^rs^ = %^ E \«sWls) ■ rs',sfS[els - e^,] , (12) 

where the 5— function is understood to be broadened to a width much larger than the 
level spacing. Assuming that there is no correlation between angular momentum matrix 
elements and orbital energy differences, it follows from the sum rule mentioned above that 

The absence of strong correlations between energy differences and angular momentum 
matrix elements, a property that we find somewhat surprising, has been verified numerically 
as we discuss below. The character of the nanoparticle quasiparticle energy spectrum changes 
when these intensive energy scales become comparable to the nanoparticle level spacing 6. 
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FIG. 1: Single-particle energy shifts caused by spin-orbit interactions in a ferromagnetic Cobalt 
nanoparticle of 143 atoms. The single-particle mean level spacing at the Fermi energy is 5 ~ 4.3 
meV. (a) Magnetization in the z-direction. (b) Magnetization in the x-direction. 

B. Numerical Results for a Co nanoparticle 

The qualitative considerations of the previous section provide a framework for thinking 
about the effects of spin-orbit interactions. Our microscopic model, on the other hand, 
allows us to explore realistic magnetic nanoparticle systems in great detail. We have studied 
numerically nanoparticles containing up to 260 atoms. Most of the results presented below 
are for hemispherical 143-atom nanoparticles. The calculations have been performed with 
a R12000 300MHz processor on an SGI Origin 2000 computer. Diagonalizations rely on 
LAPACK drivers. A single diagonalization of the Hamiltonian for a 143-atom cluster has a 
running time of approximately 1 hour, and requires around 750Mb of internal memory. In 
Fig. |]l| we plot the energy shifts caused by spin-orbit interaction for a hemispherical Cobalt 
nanoparticle of 143 atoms with a f.c.c. crystal structure. For this nanoparticle size, the 
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FIG. 2: Variation of a few quasiparticle energies of a 143-atom Co nanoparticle as a function of 
the direction of the magnetization fl. The Fermi level is the dotted line, (a) lies in the xy-plane 
and <I> is the angle with x-axis. (b) 0, lies in the zx-plane and is the angle with the z-axis. 

minority and majority single-particle mean-level spacings at the Fermi level are Si ~ 4.9 
meV and 6-^ ~ 50 meV respectively, when spin-orbit interactions are absent. The single- 
particle mean-level spacing averaged over all states (i.e. without distinguishing between 



majority and minority levels) is 5 ~ 4.3 meV at the Fermi level|^2|. The spin-orbit induced 
shifts are both positive and negative and their absolute values go from 1 meV up to 10 meV. 
The average absolute value of the energy shifts is 2.7 meV when the magnetization is in 
the z-direction, and 2 meV when the magnetization is in the xy-plane, consistent with the 
rough estimates above. We also note that there is not a strong correlation between the sign 
of the energy shift and the energy of the orbital. 

Because of the spin-orbit interaction, each individual eigenlevel has an energy dependence 
on the spin-splitting field (or magnetization) direction Cl. To illustrate typical properties of 
these dependences, we plot in Fig. 0] the variation of a few energy levels around the Fermi 
level for cases in which the magnetization rotates in 2;x-plane and a;y-planes respectively. In 
the absence of spin-orbit interactions there would be no dependence of any of these orbital 
energies on magnetization orientation. 

Notice that the angle dependence in the xy-p\ane is considerably weaker than in the 
zx-plane, for the hemispherical nanoparticle we consider. This trend indicates that for our 
nanoparticles, it is the overall shape which dominates the spin-orbit induced anisotropy 
physics. For this size cluster, there are many narrowly avoided level crossings, a point to 
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FIG. 3: Single-particle level anisotropies in the xy- and zx-planes for a 143-atom nanoparticle. 

which we return below. The difference in eigenvalue sensitivity to magnetization rotations in 
the two planes is clearly visible in Fig. 0, where we plot the single-particle level anisotropies 
^n,s{h'lz) — en,s{h%x) and en^sih^x) — ^n,s{h^y) versus the eigenvalue index n. Here and 
are the magnitudes of the self-consistent spin-splitting field when its direction is along x 
and z respectively. We note that typical change in orbital energy between x and z direction 
magnetizations is only ~ 30% smaller than the typical total shifts induced by spin-orbit 
interactions. On the other hand the typical difference in orbital energy between x and 
y direction magnetization is five times smaller than the corresponding spin-orbit induced 
energy shift. All orbital energies are relatively insensitive to the magnetization orientation 
within the xy-plane. In both cases the correlation between position within the band and the 
sign and magnitude of the energy shift is weak. In addition, energy shifts at nearby energies 
are weakly correlated. That is, the correlation function 

where the average (....) is over the occupied levels n, drops to zero very rapidly with k, as 
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FIG. 4: Correlation function of single-particle anisotropics as defined in Eq. |14]. The correlation 
function drops immediately to zero for k > 0. Note that the A; = value of the correlation function 
is equal to the width of the anisotropy distribution, plotted in Fig. . 

clearly shown in Fig. [Q. 

It is useful to consider the distribution of the quasiparticle anisotropies, P{Sen,s)- As 
an example, in Fig. ^ we plot the distribution of anisotropies in the zx-plane, 5e„^s = 
[^n,s{h*z) — en,s{h*x)], constructed with the = 9 x Afa = 9 x 143 occupied single-particle 
states of a 143-atom nanoparticle. The distribution has a width - characterized by the 
root mean square - A(5e„,s) ~ 2.2,^^/11^ ~ 2.9 meV and a much smaller mean value 
{S^n,s) = {^n,s{h*z) — en,s{h*x)) ~ IS/icV. Notc that A{6en,s) is exactly equal to the k = 
value of the correlation function displayed in Fig. 0]. Single particle anisotropies of groups 
of orbitals over within a specified energy range tend to be anti-correlated, leading to typical 
averages smaller than the variance of the distribution. The large difference between the 
distribution mean and variance will play an important role in Sec. [IV B| , where we discuss 
fluctuations of the anisotropy energy as a function of electron number. 
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FIG. 5: Distribution function of single-particle anisotropics in the zx-plane, 5en,s = ['^n,s{h* z) — 
^n,s{h*x)]^ for a 143-atom nanoparticle. The mean value is {5en,s) = 15/ieV, the width of the 
distribution (enclosed by the vertical dashed lines) is A(Je„^s) = 2.9 meV. 



The width of the distribution, A(5e„,s), gives a measure of the average magnitude of 
the single-particle level anisotropy, and the ratio A{5en,s)/^ characterizes the strength of 
mixing between quasiparticle orbitals that results from spin-orbit interactions. As mentioned 
previously, this identification of weak and strong spin-orbit interaction regimes is equivalent 
to the usual onepi], |53|, ^ based on the comparison of the spin-orbit scattering time 



Tgo and the mean- level spacing 5. 

When 5T^o/h >> 1, a limit achieved for small enough particle size for any value of ^d, 
spin-orbit coupling is a relatively weak effect and there is little mixing between spin-up and 
spin-down states. As a consequence, the level crossings between states of predominantly 
opposite spins that occur as a function of the magnitude and orientation of will be 
only weakly avoided. With increasing particle size, 5 decreases and we enter the regime 
of strong-spin orbit interaction and strong level repulsion. The single-particle spectrum 
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FIG. 6: Distribution function of spin-orbit matrix elements for a spin-polarized 143-atom nanopar- 
ticle. \n,s) and \n',s') are single-particle levels in absence of spin-orbit interaction, (a) \n,s) and 
\n',s) have like-spins; (b) \n,s) and \n,—s) have opposite-spins. Solid and dotted lines are for the 
cases of nearby levels (|e^_^ ~ ^n' s'l/*^ < 3) and any pair of levels respectively. 

becomes relatively rigid, level crossing will be strongly avoided, and 6 will limit the variation 
of individual levels as a function of the magnetization direction. 

Within our model we have found that the cross-over between these two regimes is very 
broad, it starts for nanoparticles containing of order 200 atoms and, as we argue below, it 
will be completed when the nanoparticles contain approximately 1000 atom. For a 143-atom 
nanoparticle, we are already approaching the cross-over regime; we find A{6en,s) ~ 2.9meV 
while the single-particle mean-level spacing at the Fermi level is 5 ~ 4.3 meV. In this regime 
level crossings will be moderately avoided. 

The typical size of avoided crossing gaps between opposite spin orbitals in small nanopar- 
ticles can be understood by the following argument. We first note that the unperturbed 
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orbitals satisfy the following sum rule 

^|(n',i |L -51^,1)12 = (n|L_L+|n))/4~ 2/3 (15) 

n' 

The estimate for the right hand side of Eq. [l^ is based on the on the same considerations 
leading to Eq. |Tl|, and uses that ~ 2/3L ■ L ~ 8/3. If angular momentum matrix ele- 

ments between orbitals are not correlated with the energy differences between these orbitals, 
and the matrix elements are reasonably narrowly distributed, this sum rule can be used to 
estimate the typical matrix element. Expectation values of the angular momentum opera- 
tors are zero because of angular momentum quenching, and a finite fraction of the matrix 
elements vanish because of symmetries present in our rather regularly shaped nanoparticles. 
Aside from these features, we find numerically that correlations between matrix elements 
and energy differences are too small to be clearly observable. Fig. |^ shows the distribution 
function we have obtained for matrix elements between opposite and like-spin states. We 
have considered both the matrix element distribution for closely spaced levels and the dis- 
tribution for any pair of levels, not necessarily nearby. The distributions are found to be 
very similar. Approximately 50% and 70% of matrix elements are zero for opposite-spin and 
like-spin cases respectively. 

Based on these numerical results and the sum rule Eq. we estimate the typical value 
of \{n' , I \Ha,Q\n, DP as divided by half the total number of s, p and d orbitals, 9Afa/2. 



This implies a typical non-zero matrix element equal to ~ ^(^y .W/A/l. 

In Fig. 1^ we plot the average square matrix element vs. energy difference for a 143 atom 
cluster, obtaining remarkably precise agreement with this estimate provided that the energy 
differences are much smaller than the (i-band width. 

The average square matrix element for like spin orbitals is approximately a factor of two 
smaller, consistent with the type of argument presented above which would imply propor- 
tionality to {n\Ll\n) in that case. 

We expect the dependence of quasiparticle energies on the magnitude and orientation 
of the order parameter, and also on external fields discussed below, to change in character 
when the typical matrix element becomes comparable to the level spacing, i. e. when 
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FIG. 7: Average squared matrix elements of the spin-orbit interaction vs. energy difference for 
a polarized 143-atom nanoparticle. \n,s) and \n',s') are quasi-particle states without spin-orbit 
interactions. The four curves correspond to the four possible spin combinations. 

Our numerical calculations of quasiparticle spectra are consistent with using the condition 
^'^so ~ £50 = as a criterion for the start of the crossover to the strong coupling limit, and 
the condition that the typical avoided crossing gap estimate be equal to 6, as a criterion for 
completion of the crossover to strong-coupling relatively rigid spectra. 

For Cobalt nanoparticles the later condition is reached for N'a ~ 2000; for smaller 
nanoparticles, the quasiparticles generally have somewhat well-defined spin character, some 
Poisson character in their spectral statistics, and complicated evolution patterns with exter- 
nal field and order parameter variations. For larger nanoparticles, which we are not however 
able to study numerically, we expect that quasiparticles will have strongly mixed spins, and 
more rigid spectra with smoother evolution patterns. All the ferromagnetic nanoparticles 
that we are able to study here, and many nanoparticles studied experimentally, are in the 
crossover regime. Note that since 5 oc Wd/Ma, i^d/i^-^a) oc hr^Q] the two conditions differ 
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FIG. 8: Magnetic anisotropy energies for a 143-atom hemispherical nanoparticle. The dotted hne 
represents [E{h'^ z) — E{h^ Cl)] vs. the angle Q between Q and z, when is in the zx-plane. The 
black line represents [E{h^ x) — E{h*^ Cl)] vs. the angle $ between Cl and x, when Vl lies in the 
xy-plane. 

quantitatively not parametrically. 

The total ground-state energy of the particle E{h*^tl), obtained by summing over the 
lowest N occupied orbitals, depends on the direction of the self- consistent spin-splitting 
field Vt. In Fig. we plot E{h*^z) — E{h^Q) vs. the angle B between Q and z (when 
Q lies in the zx-plane) and E{h%x) — Eih^^tt) vs. the angle $ between and x (when 
f2 lies in the 2;x-plane). From the figure we can see that the xy-plane is almost an easy 
plane for the model nanoparticle, except for a weak energy dependence which generates four 
easy axis in the directions (±x ± y)/2. Again this property reflects the large significance 
of the overall sample shape. The four easy axis directions are remnants of the magnetic 
anisotropy symmetry in bulk f.c.c. ferromagnets. Group theory considerations demand that 
a bulk f.c.c. ferromagnet have an easy axis in one of the directions perpendicular to the 
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8 (±1,±1,±1)P] planes. In a hemispherical nanoparticle the f.c.c. symmetry is partially 
lifted and the magnetization is forced to lie in the xy-plane. We can define the following 
two anisotropy-energy-per-atom constants 

k{z, x) = ^ " ' — ^ O.lSmeV , (17) 

, , Eih%x) - Eihly) , , 

fc(x, y) = \r ' ~ O.OlmeV . (18) 



As expected on the basis of the qualitative considerations of Section pi A| , the anisotropy 
per occupied orbital is much smaller than the average single-particle level shift eso (approx- 
imately 200 times smaller in both cases) due to cancellations between positive and negative 
shifts mentioned above. 

We notice that k{z,x) for our nanoparticle is /anjfer than the bulk anisotropy per atom 
^buik = 60/ieV, as one would expect because of the hemispherical shape of the sample. 
However k{x,y) is actually sma/Zer than fcbuik- This comparison should be regarded with 
caution, since it is known that accurate theoretical estimates of the magnetic anisotropy for 
bulk crystals are very delicate and agreement with experiment values even in the bulk is 
still not completely satisfactory p5| ; 

we have not evaluated the anisotropy energy that results from the bulk limit of our 
nanoparticle model and it may well not agree with experiment. Nonetheless, the small value 
that we find for the anisotropy in the xy-plane might be connected to the puzzling finding, in 
both tunneling|0] and switching- field experiments on single ferromagnetic nanoparticles, 
of anisotropy energies per atom a factor of five smaller than bulk values. 



IV. FLUCTUATIONS OF MAGNETIC ANISOTROPY 

So far we have considered only average values of the anisotropy constants. In ferromag- 
netic nanoparticles, however, these quantities are also characterized by large fluctuations as 
a function of experimentally relevant parameters. Anisotropy fluctuations are the topic of 
this section. 
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FIG. 9: Change of quasiparticle energies with increasing spin-splitting field \h\ near /i* for a 143- 
atom nanoparticle. 200 levels are shown. The Fermi level is the thick black line. The vertical white 
dashed lines indicate the positions of the charge redistributions displayed in Fig. |1C]. The vertical 
black dashed line corresponds to the self-consistent spin-splitting field. 

A. Charge- induced fluctuations as a function of the spin-splitting field 

When spin-orbit interactions are included, all quasi-particle eigenstates have mixed 
majority-spin and minority-spin character. 

Our specific calculations, however, are for small enough particles such that most eigen- 
states still have predominantly one spin character. It is convenient to use this "predominant 
spin" to label the states when discussing their dependence on \h\ or on an external magnetic 
field. In a paramagnetic system (h* = 0) in absence of external field, there is a Kramer 
degeneracy that pairs up eigenstates with opposite predominant spin character. The degen- 
eracy is lifted in the ferromagnetic state. Majority-spin states will move down in energy while 
minority-spin states will move up as \h\ increases. Because of spin-orbit-induced level repul- 
sion, all of the level crossings are avoided; the levels evolve continuously with gradually 
changing their spin and orbital character. In Fig. we plot the variation of two hundred 
levels as a function of (^s/UbI^I near QsP-Bh* = 2.2 eV, for a 143-atom nanoparticle. The Fermi 
level is the thick black line lying in a region of predominantly minority-spin quasiparticles 
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FIG. 10: Relation between charge fluctuations and anisotropy in 143-atom nanoparticle. (a) Total 
charge fluctuation defined in Eq. Charge fluctuations of the order of one electron result in 

fluctuations in the magnetic anisotropy. (b) Anisotropy energy in the zx-plane vs. the magnitude 
of the spin-splitting field taken as a free parameter. 

(lines with positive slope). By increasing \h\, individual majority-spin quasiparticle energies 
(negative slopes) come down from regions above the Fermi level, creating avoided crossing 
gaps as they approach minority-spin quasiparticles moving in the opposite direction. Since 
a 143-atom nanoparticle is still toward the weak spin-orbit coupling limit, the level crossings 
will be avoided only weakly. Whenever one of the majority-spin quasiparticles crosses the 
Fermi level, there is a change in the spin character of one of the quasiparticle levels and 
the total spin of the nanoparticle increases approximately by one. Due to the exchange 
interaction, such spin-flips bring about charge redistribution inside the nanoparticle, which 
in turns gives rise to fluctuations in the anisotropy energy. In Fig. PD| a] we plot the total 
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— * 

atomic site charge redistribution, when h = \h\ ^ h + Ah, for Ah = 12.5 meV 

spih) = Y.\p,{h + Ah)-piih)\^Yl 

i i 

where pi{h) is the total charge at atom i. It is seen that 6p{h) changes by one in the small 
interval Ah = 12.5 meV where majority-spin quasiparticles weakly avoid crossing the Fermi 
level. (Note that this is a charge redistribution, not a change in total charge.) This result 
is not unexpected, since the majority and minority spin orbitals should have uncorrelated 
spatial distributions. The charge redistribution is likely overstated by our simplified model, 
however, since it does not account for long-range Coulomb interactions and related screening 
effects. It is, however still a relatively small ~ 1/A^ effect. The anisotropy energy fluctuations 
associated with the level crossings that occur as a function of h, shown in Fig. |IO| b] are much 
larger in relative terms. 

Variation in h corresponds to variation in the amplitude of the magnetic order parameter. 
In micromagnetic modeling of nanoparticles, it is implicitly assumed that the magnitude of 
the order parameter is flxed and that this collective degree of freedom can be ignored in 
modeling nanoparticle properties. The substantial dependence of anisotropy energy on h 
that we flnd demonstrates that the amplitude and orientation fluctuations of the order 
parameter can be strongly coupled in small magnetic nanoparticles. Related anisotropy 
fluctuations as a function of electron number are important in understanding the addition- 
potential-spectroscopy single-electron-transistor experiments]^] which are the topic of the 
next section. 



dpi{h) 



dh 



Ah 



(19) 



B. Mesoscopic fluctuations as a function of electron number 

The analysis of the field-dependence of tunneling resonances in experiments on magnetic 
nanoparticles suggests that the anisotropy energy fiuctuates significantly from eigenstate to 
eigenstate|0, |^. In Refs. 0, |^ the effects of these fluctuations were mimicked by using two 
different anisotropy energy constants and kN±i = + 6k± for N- and ± 1-electron 
states. By assuming 6k±/kN in the range of a few percent, it was possible to explain the 
non-monotonic behavior of the tunneling resonances seen experimentally. 

Within our microscopic model, we have calculated k^ and k]\f±i as deflned in Eqs. (^^8|). 



In light of results of the type summarized in Fig. 0, it suffices to compute the total energy 
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for two different directions of the spin-splitting field and take the difference. The total energy 
for a A^-electron system, for spin splitting fields along directions Cli , i = 1, 2, is given by: 



N 



{h: 



★ \2 



(20) 



n=l 



where h* is calculated self-consistently for a given fixed direction. The anisotropy-energy 
per-atom constant is then 



(21) 



where the subscript N emphasizes the fact that the constant refers to a A'"-electron system. 
It turns out that the value of h* depends very weakly on Q. For example, in a 143-atom 
nanoparticle \h* — h^l/h*^ ~ 10~^. This property reflects the large ratio between the total 
magnetic condensation energy and the anisotropy energy, even for nanoparticlcs. The strong 
coupling between amplitude and orientation fluctuations mentioned above, occurs only when 
a small change in h leads to a crossing between majority and minority spin orbitals. Evalu- 
ating E{hQ,^ ai h — h\ — + Sh*, expanding in powers of Sh*/h\ and remembering that 
/i2 minimizes E{hri^ yields 



n=l 



N ^6h\2 



-I N r 7W- 



(22) 



for any N. 

Adding an electron to the system changes the magnitude of the spin-splitting held h* 
h*', but again \h*' — h*\/h* « 1. Thus we obtain 



+ 



n=l 



N + 1 f5h\2 



N + l ^Sh*Y 



Replacing 



€N+i{h*^ ^i) - eN+i{hl O2) A[(5ejv+i(f^i, ^^2)] 



(23) 



(24) 



where A[(^ejv+i(f^i, ^^2)] is the width of the distribution of the single-particle anisotropics in 
plane containing the two directions Cli and Q2, we flnally obtain 



(25) 



25 



z 0.1 




(gUglhl) [eV] 




(gl^iglhl) [eV] 



FIG. 11: Relative shift 5k±/kN = (^7V+i — kj\f)/k]\f of the anisotropy energy constant k]\f{z,x), 
when the electron number increases by one, ^ + 1, as a function of spin-splitting field, 
(a) 260-atom nanoparticle. (b) 143-atom nanoparticle. The self-consistent spin-splitting field is 
h* ~ 2.2 eV for both a 143-atom and a 260-atom nanoparticle. 

The fluctuations of fc^v due to an additional electron are therefore of the order of 
A[SeN+i{^i,^2)]/J^a, and are regulated by the mean level spacing, 6, which suppresses 
the magnitude of A[6eN +1(^1,^2)] at large nanoparticle sizes as we discuss below. 

For a 143-atom dot, A[(5eAr+i(i7i, i72)] ~ 2.9meV in the zx-plane, where ki^{z,x) ~ 0.13 
meV. Therefore 5k±/k^{z,x) = -^A[6eM+i{^i,^2)]/kN{z,x) ~ 15%. These estimates are 
confirmed by a direct calculation of kj\f{z,x) and ki\i^i{z,x) for 143-atom and 260-atom 



nanoparticles, shown in Fig. |[TT[]. Here kN{z,x) and kN+i{z,x) are plotted as a function of 
QsfJ'Bh, taken as a free parameter (the self-consistent value is close to gsfJ^sh* = 2.2 eV for 
both a 143-atom and a 260-atom nanoparticle). For a 143-atom nanoparticle the calculated 
anisotropy energy fluctuations are of the order of 20%. For a 260-atom nanoparticle the 
anisotropy energy fluctuations are smaller, but still close to 5%. 
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We can understand the surprisingly large fluctuations in anisotropy energy with parti- 
cle and electron number by the following observations. Our numerical results (see Fig. |^ 
and Fig. |^) indicate that the contribution of a given orbital to the anisotropy energy is 
chosen essentially at random from a distribution which, for the zx-plane anisotropy, has 
a width A{6en,s) ~ '^■'^ ^df^d ~ 2.9 meV and a much smaller mean value {Sen,s) = 
{en,s{h*z) — en,s{h*x)) = AfakN{z,x)/N ~ 15/ieV, where kN{z,x) ~ 0.13 meV. The rela- 
tive change in the anisotropy energy expected when one additional orbital is occupied is 
~ 4:.3^yWd/{kN{z,x)Afa), which can easily be larger than 1% for Afa < 1500. 

We expect strong level repulsion in the rigid spectrum of larger nanoparticles to be 
accompanied by more regular behavior of the anisotropy energy per atom, with smaller 
variations as a function h*, and electron number. 

We note that fluctuations in the contribution of a given orbital to the anisotropy energy 
should not exceed ~ 6, which vanishes in the limit of very large particles. We expect that the 
distribution function of contributions to anisotropy from orbitals near the Fermi energy to 
become narrow when 6 is smaller than the average anisotropy energy contribution J2 lO/xeV. 
This condition is satisfied for particles containing more than ~ 10^ atoms. 



C. Anisotropy energy dependence on particle atom number and shape 

We conclude this section with a few remarks on the dependence of the magnetic anisotropy 
energy on nanoparticle size and shape. In Fig. we plot the anisotropy energy per 



atom for nanoparticles of different sizes. For a small number of atoms (below 60) the 
anisotropy per atom is very large and decreases rapidly with size. In this regime, non- 
extensive surface contributions, which are present because of the abrupt truncation of the 
lattice, are clearly playing a dominating role. By increasing the size of the nanoparticle, we 
must eventually reach a regime where total magnetic anisotropy becomes proportional to 



the nanoparticle volume. Fig. [|T^] shows that this regime is not yet fully established even for 
nanoparticles with hundreds of atoms, and fluctuations are still pronounced. In this regime 
the anisotropy per atom is still 2-3 times larger than the bulk value for Cobalt, which is 
0.06 meV. There is no regime in which bulk (proportional to ACi) and surface (proportional 
to N'a^^) contributions to the anisotropy can be cleanly separated. 

In is clear that for small particles the shape of the nanoparticle plays an important role 
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FIG. 12: Anisotropy energy as a function of nanoparticle size for two values of the spin-orbit 
coupling (^ci- A splitting field of 2 eV is assumed for all sizes. 



in the determination of the magnetic anisotropy. For example, as seen in Fig. ||13[, for 
spherical particles we find anisotropy energies which are one order of magnitude smaller 
than the zx anisotropy of hemispherical particles, comparable instead to the xy anisotropy 
of hemispherical particles with the same number of atoms. A 140-atom spherical particle 
has an anisotropy energy per atom ^ 0.01 — 0.02meV. 

It is interesting to notice that both tunneling]^ and switching- field experiments on 
single ferromagnetic nanoparticles find anisotropy energies per atom which are of the order 
O.OlmeV, a factor of five smaller than the bulk value. The nanoparticles in tunneling 
experiments are roughly hemispherical |]T], |^, whereas the nanoparticle shape in Ref. |0] 
is close to spherical. 

Further theoretical studies that focus on the relationship between nanoparticle shape and 
anisotropy could be helpful to efforts to engineer ferromagnetic nanoparticles whose shape, 
size, and crystal structure are tuned to produce desired magnetic properties. 
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FIG. 13: Magnetic anisotropy energy in the zx-plane, [E{h'^ z) — E{h^ Q)] for spherical nanopar- 
ticles of two sizes. is the angle between z and Cl. By symmetry, the anisotropy energy in the 
xy-plane is now the same as in the zx-plane. The spin-splitting field is assumed to be 2 eV for 
both sizes. 

V. HYSTERESIS AND VARIATION OF SINGLE-PARTICLE LEVELS IN AN EX- 
TERNAL MAGNETIC FIELD 



In the last section we shall investigate the effect of an external magnetic field on the total 
energy of the nanoparticle. This will allow us to make a connection between our microscopic 
model and more familiar classical micromagnetic energy functional expressions which are 
normally used to interpret the results of switching- field experiments[0]. We shall also study 
the dependence of the quasiparticle energy levels on external magnetic fields; this complex 
behavior has been probed directly in recent tunneling experiments 0]. 
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FIG. 14: "Uniaxial" anisotropy energy as a function of the direction of the spin-sphtting field 0, 
for different magnitudes of the external magnetic field. The external field is in the zx-plane at an 
angle 9 = tt/4: with the z-axis. G is angle between Q (lying in the zx-plane) and z. 

A. Hysteresis 

Ideally one would like to study the three dimensional energy landscape E{h* , He^t) as a 
function of f2(G, $) = h*'/\h*\, for several values of -ffext- Here we consider two simplified 
cases. In Fig. we plot the total energy difference E{h\z) — E{h^{l) as a function of an 
external magnetic field i^ext- Here Cl is allowed to rotate in the 2;a;-plane. In this picture 
we recognize the familiar features of a micromagnetic energy functional characterized by a 
uniaxial magnetic anisotropy: at ifext = there are two degenerate minima, separated by 
an energy barrier, representing two equivalent easy directions ±x. When a magnetic field 
is applied in the 2;a;-plane, the degeneracy is lifted and the minimum in the —x direction 
becomes a metastable local minimum, until a switching field is reached at which the local 
minimum disappears. This happens at |-ffext| ~ IT for a 143-atom nanoparticle. This 
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FIG. 15: Variation of the magnetization direction of the stable minimum of Fig. |14|, as a function 
of l-ffcxtl- is the angle between Cl and z. i^cxt is in the zx-plane at 7r/4 from the z- axis. 

simple model captures the essence of classical hysteresis in ferromagnetic nanoparticles. 
Quantum mechanically, for sufficiently weak fields the system can still sit in the quantum 
state characterized by a magnetic moment pointing along the classically metastable direction 
until the switching field is reached, the system relaxes to its true ground state, and the 
magnetization orientation collective coordinate changes discontinuously |^ . 

With a further increase of the external field, the magnetization is gradually twisted from 
the easy x-axis toward the direction of the magnetic field, as shown in Fig. [|15[. At low 
fields (l-ffextl < IT), the direction of the magnetization corresponding to this local minimum 
is very close to the xy-plane. 

Because our hemispherical nanoparticles have a nearly isotropic easy xy-p\a.ne, the mag- 



netization will rotate in the xy-plane in response to weak applied external fields This is 
exemplified in Fig. |]16i, where we plot the the total energy as a function of l-f^extl? when Q 



lies in the xy-plane. |-ffext| is oriented in the direction {6 = 7r/4, = 7r/16). In the absence 
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FIG. 16: Anisotropy energy in the xy-plane for different magnitudes of an external magnetic field. 
The direction of the external field is = 7r/4, = 7r/16). $ is the angle between Cl (lying in the 
xy-plane) and x. 

of the external field, the ground state of the nanoparticle is four-fold degenerate, with the 
degeneracy corresponding to one of the four magnetization directions ±x ± y, associated to 



the four local minima of Fig. As the external field increases, three of these minima 

will become classically metastable with small barriers separating them from each other and 
the true ground state. In general the different local minima will lose their metastability at 
different external field strengths. In the case we consider, the first switching field is reached 
at |-ffcxt| ~ 0.15T, when the minimum originally at —x — y disappears. If the system starts 
out in this local minimum at zero external field, at the switching field it will jump to the 
ground state with a corresponding rotation of the magnetization in the x + y-direction. Since 
everything takes place essentially within the xy-plane, where the typical anisotropy energies 
are one order of magnitude smaller than in the 2;x-plane, the scale of the coercivity is also 
much smaller than the one in the zx-plane. More complicated hysteretic behaviors are also 
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possible, in which, by starting from a different minimum at zero external field, the system 
jumps first from one metastable to another metastable state and only at a second switching 
field reaches its true ground state. 



B. Dependence of quasiparticle levels on external fields 

The hysteretic behavior that we have seen in the ground-state properties of a ferro- 
magnetic nanoparticle has profound implications for the magnetic field dependence of its 
low-energy elementary excitations. In ferromagnetic metals there exist two kinds of ele- 
mentary excitations: collective spin excitations associated with magnetization orientation 
degrees of freedom, and particle-hole excitations. In a ferromagnetic nanoparticle the dis- 
tinction between these two kinds of excitations is partly obscured by the effect of spin- 



orbit coupling [|TT|, p6| . Sorting out this problem is particularly important in order to inter- 
pret current tunneling experiments in single-electron transistors]^. Here we examine only 
particle-hole excitations around the Fermi level, which are immediately available within our 



Hartree-Fock treatment of the many-body Hamiltonian. In Fig. |]T^ we plot the magnetic 
field dependence of a few single-particle energy levels around the Fermi energy. As a simpli- 
fied illustration, the levels are calculated assuming that the ground-state dependence on the 



external field is as described in Fig. 0]: at low field the ground-state magnetization is ori- 
ented around the — x-axis until the switching field is reached, whereupon the magnetization 
is reversed along a direction around the x-axis. The corresponding field dependence of the 
quasiparticle is quite complex. In the small- field regime (|-ffext| < 1 T) there is a hysteretic 
switching at a certain |-ffext| = Hsw ~ 1 T, due to an abrupt change of the ground-state 
magnetic moment. There is basically no correlation between single particle states before and 



after reversal. Notice in particular that the levels can jump either up or down at ifswPS 
Furthermore the quasiparticle energies have continuous non-monotonic variations, which 
seem to differ randomly from level to level. In the large field regime (|i/cxt| >> H^^), the 
quasiparticle energies depend roughly linearly on |-ffext| and their slopes almost all have the 
same sign. In the small-field behavior, the variation of the quasiparticle energies as a func- 
tion of the external field is due to a combination of the rotation of the particle's magnetic 
moment direction Q, and Zeeman coupling. We have shown above (see Fig. |^) that the 
dependence of the quasiparticle energies on Q varies randomly from level to level. Thus the 
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FIG. 17: Variation of quasiparticle energy levels e„^s(/i*, //ext) around the Fermi level in an external 
magnetic field. The levels are calculated assuming that the ground-state of the nanoparticle changes 
in the external field in the way described in Fig. |14]. 



complex non-monotonic behavior at small fields can be understood within our model. At 
large fields (|-ffext| >> -f^sw) ^ is slowly twisted toward the field direction and the Zeeman 
coupling plays the dominating role in the field dependence. Moreover, because the Fermi 
level lies in a region of predominantly minority-spin energy levels, it is expected that almost 
all particle and hole excitations around the Fermi level have the same high-field slope 0]. 

In the discussion above we have assumed, out of simplicity, that the magnetization is 
arbitrarily constrained to rotate in the ^x-plane. At weak external fields, however, the 
magnetization will stay close to the xy-plane. For a weak external field in the {6 = 7c/4,(f) = 
7r/16) direction, the actual dependence of the ground-state energy on the magnetization 
direction at different field strengths will be the one represented in Fig. |T6[. Consequently, 
the change of quasi-particle energies as a function of the external field will reflect this more 
complicated behavior. We illustrate this point in Fig. |l8|-a] where, as an example, we 
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FIG. 18: Variation of quasiparticle energy levels en^sih* , H^xt) around the Fermi level in an external 
magnetic field, (a) The levels are calculated assuming that the ground-state of the nanoparticle 
changes in the external field in the way described in Fig. ||l6|. (b) For magnetic fields larger than 1 
T, the magnetization starts to develop a non-negligible component perpendicular to the xy plane, 
while its projection in the xy-plane is parallel to the component of -ffext in that plane. 

assume that at zero field the magnetization points in the —x + y direction. As the magnetic 
field increases, there will be a first jump in the in the quasiparticle-level dependence at 
^ 0.15 T, when the magnetization re-orients itself in the (-l-a; -|- y)-direction. There will 
a second hysteresis jump in the quasi-particle energies at ifg^ ^ 0.3 T, when the magnetic 
moment finally switches to the (-l-a; — y)-direction, corresponding to the only metastable 
configuration left at this field. For a different starting local minimum at |-ffext| = 0, the 
variation of the quasiparticle energies as a function of the external field will be in general 
different. By increasing the magnetic field beyond 1 T, the magnetization will start to 
develop a non negligible component perpendicular to the xy-plane, while its component in 
the xy-plane will become essentially frozen along the direction of the component of i^ext in 
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that plane. 

The quasiparticle-energy dependence on an external magnetic field presented here, has 
striking similarities with the field dependence of the tunneling resonance energies in single- 
electron tunneling experiments [jT], |^. Notice that for the hysteretic behavior occurring as a 
result of the rotation of the magnetization in the vicinity of the easy xy-plane, the order 
of magnitude of the coercivity that we find is close to the experimental value. The only 
discrepancy between the results of our theoretical model and the experiment is the mean 
level spacing 6^es of the low-energy excitations. Since in our model we are only taking into 
account Hartree-Fock quasiparticles, the expected level spacing of the low-lying excitation 
is approximately 6i, which is a few meV in our 143- atom nanoparticle and should be of the 
order of 0.5 meV for a 1500-atom nanoparticle considered in the experiment. 

This value is still larger than the value 5res < 0.2 meV observed experimentally. We believe 
that a unified and consistent inclusion of collective spin excitations, and possibly also the 



non-equilibrium transport effects proposed in Refs. M, 10 1, could resolve this confusion. 



VI. CONCLUSIONS 



In this article we have investigated the effects of spin-orbit interactions on the proper- 
ties of ferromagnetic metal nanoparticles. In particular, we have focused on their novel 
microscopic magnetocrystalline anisotropy physics, and on hysteresis in the quasiparticle 
excitations spectra of metallic nanomagnets. Our analysis, based on qualitative consid- 
erations backed up by numerical studies of a generic model for ferromagnetic transition 
metal nanoparticles, provides an understanding of some emergent properties of their quasi- 
particle states. We find two regimes separated by a broad crossover and characterized by 
the comparison of several characteristic energy scales. For small nanoparticles with fewer 
than A/'a ~ 200 atoms, the single particle mean-level spacing 6 is larger than spin-orbit in- 
duced energy-shifts eso in the quasi-particle spectra. These shifts have the same typical size 
as the spin-orbit scattering lifetime broadening energies of very large particles, hr^Q. The 
quasi-particle levels of small nanoparticles in which St^o > ^ evolve in a complicated way 
as a function of magnetization orientation and external magnetic field, with relatively small 
avoided crossing gaps and spin-orbit shifts of nearby orbitals that are nearly uncorrelated. 
The size of the avoided crossing gaps is determined by matrix element of the spin-orbit 
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coupling operator between quasiparticle energy levels that are adjacent in energy. Surpris- 
ingly, even though expectation values of these matrix elements vanish because of angular 
momentum quenching, typical values between energetically adjacent orbitals are compara- 
ble to those for orbitals at any position in the spectrum. Eventually, for nanoparticles with 
more than ~ 1000 atoms, the typical avoided crossing gap estimated in this way becomes 
comparable to the level spacing and we expect the crossover to the strong coupling limit is 
complete. Nanoparticles with fewer than Ma ~ 1000 atoms can easily be prepared with cur- 
rent synthesis techniques and systems of experimental interest are often in the middle of the 
crossover between small particle and bulk (weak and strong spin-orbit coupling) limits. For 
example, the nanoparticles investigated by electron tunneling experiments contain between 
50 and 1500 atoms 0, 0]. 

For nanoparticles in the size range we are able to study numerically. Ma smaller than a 
few hundred, we find that the anisotropy energy per atom displays large changes of order 
several percent when the electron or atom number changes by one. Our analysis allows 
us to make a connection between the microscopic model of a metal nanomagnet and more 
familiar classical micromagnetic energy functional expressions which are normally used to 
interpret the results of switching-field experiments^. The ground-state energy as a function 
of the magnetization direction is characterized by minima separated by energy barriers. The 
quasi-particle levels exhibit a complex non-monotonic behavior and abrupt jumps when the 
magnetization direction is reversed by an external magnetic field. The nanoparticles inves- 
tigated by electron tunneling experiments contain between 50 and 1500 atoms[|l], 0. The 
results that we have presented here are therefore particularly relevant for the interpretation 
of these experiments. In particular, we find that the anisotropy fluctuations inferred from 
interpretations of these experiments are indeed to be expected. Similarly, the dependence of 
the tunneling resonances on the external magnetic field is qualitatively similar to the behav- 
ior of the quasi-particle excitations of our model, although it appears necessary to invoke 



non-equilibrium quasiparticle configurations |]T0| in the nanoparticle in order to understand 
the size of the quasiparticle level spacing. 
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